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Abstract 

Quantum mechanics of a general one dimensional dissipative system 
investigated by it's coupling to a Klein-Gordon field as the environ- 
ment using a minimal coupling method. Heisenberg equation for such 
a dissipative system containing a dissipative term proportional to ve- 
locity obtained. As an example, quantum dynamics of a damped 
harmonic oscillator as the prototype of some important one dimen- 
sional dissipative models investigated consistently. Some transition 
probabilities indicating the way energy flows between the subsystems 
obtained. 



1 Introduction 

In classical mechanics dissipation can be taken into account by introducing a 
velocity dependent damping term into the equation of motion. Such an ap- 
proach is no longer possible in quantum mechanics where a time-independent 
Hamiltonian implies energy conservation and accordingly we can not find a 
unitary time evolution operator for both states and observable quantities 
consistently. 
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To investigate the quantum mechanical description of dissipating systems, 
there are some treatments, one can consider the interaction between two sys- 
tems via an irreversible energy flow [1,2], or take a phenomenological treat- 
ment for a time dependent Hamiltonian which describes damped oscillations, 
here we can refer the interested reader to Caldirola-Kanai Hamiltonian for a 
damped harmonic oscillator [3]. 

H (t) = e- 2 ^^- + e 2 ^-mu 2 q 2 . (1) 
W 2m 2 H K ' 

There are significant difficulties about the quantum mechanical solutions of 
the Caldirola-Kanai Hamiltonian, for example quantizing such a Hamilto- 
nian violates the uncertainty relations or canonical commutation rules and 
the uncertainty relations vanish as time tends to infinity. [4,5,6,7,8] 
In 1931, Bateman [9] presented the mirror-image Hamiltonian which con- 
sists of two different oscillator, where one of them represents the main one- 
dimensional damped harmonic oscillator. Energy dissipated by the main 
oscillator completely will be absorbed by the other oscillator and thus the 
energy of the total system is conserved. Bateman Hamiltonian is given by 

with the corresponding Lagrangian 

L = mix + jr(xx — xx) — kxx, (3) 

canonical momenta for this dual system can be obtained from this Lagrangian 

as 

dL ^ (3_ dL . (3 

V = 77 = mx - -x, p = = mx + -x, (4) 

ox 2 ox 2 

dynamical variables x,p and p, x shoud satisfy the commutation relations 

[x,p] = i, [x,p]=i, (5) 

however the time-dependent uncertainty products obtained in this way, van- 
ishes as time tends to infinity. [10] 

Caldirola [3,11] developed a generalized quantum theory of a linear dissipa- 
tive system in 1941 : equation of motion of a single particle subjected to a 
generalized non conservative force Q can be written as 

d , dT dT 8V , 
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where Q r = —f3(t)J2 a rjQj, and a rj 's are some constants, changing the vari- 
able t to t*, using the following nonlinear transformation 

t* = X (t), dt = cf>(t)dt*, 0(t) = e /oW, (7) 

together with the definitions 

9* = ^, L* = L(q,q*,n P* = |p (8) 



the Lagragian equations, can be obtained from 

d ,dL\ dU 



dt* dq* dq 



0. (9) 



where H* = J2p*Q* ~ L*. Canonical commutation rule and Schrodinger 
equation in this formalism are 

[q,P m ]=i, H ** = i W> (10) 

but unfortunately uncertainty relations vanish as time goes to infinity. [10] 
Perhaps one of the effective approaches in quantum mechanics of dissipa- 
tive systems is the idea of considering an environment coupled to the main 
system and doing calculations for the total system but at last for obtaining 
observables related to the main system, the environment degrees of freedom 
must be eliminated. The interested reader is referred to the Caldeira-Legget 
model [12,13]. In this model the dissipative system is coupled with an envi- 
ronment made by a collection of N harmonic oscillators with masses m n and 
frequencies u n , the interaction term in Hamiltonian is as follows 

N N 2 

tf'—sEc^ + g 2 ^— (11) 

n=l n=l Z " L n^ n 

where q and x n denote coordinates of system and environment respectively 
and the constants Cfi clX6 called coupling constants. 

The above coupling is not suitable for dissipative systems containing a dis- 
sipative term proportional to velocity. In fact with above coupling we can 
not obtain a Heisenberg equation like q + u 2 q + (3q = £(t), for a damped 
harmonic oscillator, consistently. In this paper we generalize the Caldeira- 
Legget model to an environment with continuous degrees of freedom by a 
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coupling similar to the coupling between a charged particle and the electro- 
magnetic field known as the minimal coupling. In sections 2, the idea of a 
minimal coupling method is introduced .In section 3 the quantum dynamics 
of a damped harmonic oscillator is investigated completly and some transi- 
tion probabilities indicating the way energy flows between subsystems, are 
obtained. In section 4, quantum dynamics of the oscillator environment is 
investigated. 



2 Quantum dynamics of a one dimensional 
damped system 

Quantum mechanics of a one dimensional damped system can be investigated 
by introducing a reservoir or an environment that interacts with the system 
through a new kind of minimal coupling term. For this purpose let the 
damped system be a particle with mass m influenced by a external potential 
v(q). we take the total Hamiltonian, i.e., system plus environment like this 

H = {p-Kf + + 
2m 

where q and p are position and canonical conjugate momentum operators of 
the particle respectively and satisfy the canonical commutation rule 

[q,p]=i, (13) 
and Hb is the reservoir Hamiltonian 

/+oo 
d*koj % b\{t)b % {t), u s =\k\. (14) 

-oo 



Annihilation and creation operators b%, &t, in any instant of time, satisfy the 
following commutation relations 

[& £ (t)4(*)] = (15) 

and we will show later that reservoir is a Klein-Gordon type equation with 
a source term. Operator R have the basic role in interaction between the 
system and reservoir and is defined by 

/+oo 
d S k[f(^)bj;(t) + r(^)bl(t)l (16) 
-oo 



4 



let us call the function f(ujj:), the coupling function. It can be shown easily 
that Heisenberg equation for q and p leads to 



q = i[H, q] 



p- R 



p = i[H,p] = - 



rn 
Ov 

dq 



(17) 



where after omitting p, gives the following equation for the damped quantum 
system 



mq 



oq 



Using (JT5)l the Heisenberg equation for b^, is 

k = b n ] = -iuJ K b n + iqf*{uj u ), 
with the following formal solution 

hit) = b$S)e** + i/> S ) f dt'e-^-*H{t'l 

J 

substituting 6g(t) from (j2(Jj) into (|18|) one can obtain 

ft Qy 

mq + ^ dt'q(t')j(t -t') = -—+ £(t) 

POO 

7(£) = 8tv du^\f{u^)\ 2 u\cosu}^t 



(18) 
(19) 
(20) 



+oo 



(21) 



It is clear that the expectation value of in any eigenstate of Hb, is zero. 
For the following special choice of coupling function 



\ 



(22) 



equation (|2ip takes the form 



dv 



mq + f3q = -— + £(£) 



4tt 2 



d 3 k 



'UJl 



(%(0)e~ i ^ i -4(0)e i ^ t ), 



(23) 



In the following we investigate for example one dimensional harmonic oscil- 
lator. 
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3 quantum mechanics of one dimensional damped 
harmonic oscillator 



3.1 quantum dynamics 



For a one dimensional harmonic oscillator with mass m and frequency ui we 
have v(q) = ^muo 2 q 2 and therefore we can write (|23|) as 



-, 0., 2 ffl 

q H q + uj q = 

m m 



(24) 



with the following solution 

q(t) = e-^{Ae iuit + Be***) + M{t) 
M(t) = i 



d A k t 



%(0) 



4(o) 



uo 2 - ul + ^ 
k m K 



(25) 



where o>i 
tions 



or 



4m- 



-. Operators A and 5, are specified by initial condi- 



A + B 

zl + luJl) A + (ZE- l u l )B = g(0)-M(0) 
2m 2m 



g(o)-M(o), 

g(0) - M(0) 
p(0) - R(0) 



m 



M(0), (26) 



solving above equations and substituting A and B in (J2*5*j) one obtains 

q(t) = e~^{ ^ - sin u;ii + g(0) cosa;it H — — g(0) sinc^t 

i?(0) . /3M(Q) . w/n , M(0) . , 

y sin wit - ^ v y sin t^t - M(0) cos u^i — sm wit} + Af (t) , 



2mui 



(27) 



6 



also substituting q(t) from (|27jl in (J2U|) we can obtain a stable solution for 
in £ — >• oo as 



Hit) = bi(0)e 



k — I. 







\ 4ir 2 u 3 , (u 2 -u£- ^a;, 

^ k v k m k 



-{uj 2 q(0) + iu % 



p(0) - i?(0) 



-M(0V 2 -iu;gM(0)} 

z/3 . 3 fep(O) sin (a, g~^ f) t -«he+V>* 



k' m k' 9 



-e 2 



ftt (o) sin^^t ^- g >' 

e 2 }, 



(21 



now substituting b^{t) from in (fTfi)) and using (|T7|). one obtains p 
mg + i?. 

A vector in fock space of reservoir is a linear combination of basis vectors 



\N{k x ),N{k 2 ), ...) B = 1 hl = 1 \0) B 



'Nihy.Nihy.... 

where are eigenstates of Hg and the operators and bt act on them as 



(29) 



b % \N(h), N(k 2 ), ...N(k), ...) B = yJNifylNih), N(k 2 ), ...N(k) - 1, ...) B 

bi\N(h), N(k 2 ), ...N(k), ...} B = \fN(k) + l\N(kx),N(k 2 ), ...N(k) + 1, ...) B 

(30) 

If the state of system in t = is taken to be \if}(0)) = \Q)b <§) \n) u where |0)# 
is vacuum state of reservoir and \n) u an excited state of the Hamiltonian 
H s = £ — h ^mu 2 q 2 , then it is clear that 



{mi ^~ + ^V(0)|V>(0)> = (n + 1) W> 



(31) 
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On the other hand from (|27|) . (|28|) and we find 

^^[(^(0)1 : X -mq 2 + ^mu 2 q 2 : |^(0))] = 0, 

Zzm^KVKO)! : ^- + im^V : |V>(0)>] 

^|?^ , f, B ]V(0))n 



27T 2 m J-oo x co 2 — x 2 + z — x 

where : : denotes normal ordering operator. Now by substituting b^(t) 
from ()28p into (|14jl . it is easy to show that 

^77^4(^(0)1 : # B (*) : |V(0)>] 



vr ./o (a; 2 - x 2 ) 2 + Kx 2 " vrm 2 7o (u; 2 - x 2 ) 2 + -^rX 2 

/9cj 3 . 1 . f 00 a?x /3a; , 1 . r°° x 2 dx 

(n+-) I — — ^- + — (n+ 1 ' 



7rm v 2 7 7o (a; 2 - x 2 ) 2 + ^z 2 7rm v 2 7 7o (w 2 - x 2 ) 2 + ^x 2 ' 

(33) 

For sufficiently weak damping that is when (3 is very small , the integrands 
in (JBlIj) have singularity points x = ±(u>i ± and by using residual calculus 
we find 

Um^ooMMlM : H B : |^(0)>] = (n+ (34) 

Comparing (|3*T|) and (|34p. one can show that the total energy of oscillator has 
been transmited to the reservoir and according to (J3*2~)) . the kinetic energy of 
oscillator tends to zero. 

If the state of system in t = is p(0) = ® |S% where = — 5 — _ Hb is 

Tr B {e~F?r) 

the Maxwell-Boltzman distribution and is an arbitrary state of harmonic 
oscillator , then by using of TrB[bUo)b^(Q)p] 3 ] = & ^Z k ^ one can show the 

eKT-l 

expectation value of kinetic energy of oscillator in t — > oo tends to 

1 -2/jA , ^ 2 2/.\ \ /"°° 



lirru-^ooi'- -mq (t) H — mcj o (£) :) = / ^ <ix 

X p2 01/ _s_ TT rfx ( 35 ) 



nm 2 Jo [(u 2 - x 2 ) 2 + ^2X 2 ](e7fr - 1) 
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3.2 Transition probabilities 

We can write the Hamiltonian (|T2|) as 

H = H + H', 
H Q = (a^a+^)u + H B , 
v R 2 

H'=-*-R+" (36) 
m 2m 

where a and are annihilation and creation operators of the Harmonic 
oscillator. In interaction picture we can write 

ai (t) = e iHot a(0)e- iHot = ae~ iu \ 

hi(t) = e^\(0)e-^ = b s (Q)e-*>*, (37) 

the terms ^p and ^- are of the first order and second order of damping 
respectively, therefore, for a sufficiently weak damping, ^ is small in com- 
parison with —p. Furthermore ^ has not any role in those transition prob- 
abilities where initial and final states of harmonic oscillator are different, 
hence we can neglect the term in H'. Substituting aj and from (|37j) 
in — — p, one can obtain Hj in interaction picture as 

H>1 = /"7 d3fc (/W° t6 jf( ) ei(w " Wt)< " /*K)o6+(0)e-^ w - w ^ 

-/(u; fe )a%(0)e- i K-+-)* + f*(u k )a%(0)e i( -^+^), (38) 

the terms containing just a%(0) and a^fet(O) violate the conservation of en- 
ergy in the first order perturbation, because ab^(0) destroys an excited state 
of harmonic oscillator while at the same time destroying a reservoir excita- 
tion state and a^bt(0) creates an excited state of harmonic oscillator, while 
creating an excited reservoir state at the same time, therefore we neglect the 
terms involving ab^(0) and a^&t(0) because of energy conservation and write 
Hj as 

H'j = -*^/_ + Jd 3 M/(^ (39) 
Now the time evolution of density operator in interaction picture is [14] 

pi(t) = C/j(Mo)p/(*o)c4(Mo), (40) 
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where Uj is the time evolution operator, which in first order perturbation is 

UAt, t Q = 0) = 1 - i C dhH'Atx) = 
Jo 

I 77] r+oo i(u— o> E )t 



(LO — LJr) , 

-r^ab^e-^}^-^!. (41) 

2 

Let pi(0) = \n) u u (n\ ® \0)b b(0| where |0)# is the vacuum state of the 
reservoir and \n) u , an excited state of the harmonic oscillator, then by sub- 
stituting Uj(t, 0) from (|4*T|) in ()40|) and taking trace over reservoir parameters 
we obtain 

p sI {t) := Tr B {p I {t)) = \n) u u (n\ 



+°° , , , , ,oSin v p ' t 



oo 



(42) 

where we have used the formula Ttb^j^b B^-g\] = 5{k — k') . For very large 

sin 2 ("p-") f 

time, we can write — zr^S — = 2irt5(uj^ — uj) which leads to the following 

( 2 ) 

relation for density matrix 

4:7r 2 u 3 nt\ f(uj)\ 2 , v / \ 

Psi{t) = \n) u M + J 1 Jl \n - 1). u (n - 1|, (43) 

from density matrix we can calculate the probability of transition \n) u — > 
\n — 1)^ as 

r n ->n-i = TV[(|n - 1)^ u (n - l\)p{t)] = 

mm v in M 4:7T 2 UJ 3 nt\ f (uj)\ 2 , . 

Tr s [(\n - l) w w (n - l|)p./(t)] = [ ^ L , (44) 

where Tr s denotes taking trace over harmonic oscillator eigenstates. For the 
special choice (|22j) . above transition probability becomes 

= 45 

m 
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Now consider the case where the reservoir is an excited state in t = for 
example p 7 (0) = \n) u u (n\ <g> ...1^)b -lft | where ...lft) B de- 

notes a state of reservoir that contains j quanta with corresponding momenta 
Pi, ...Pj, then by making use of 

rr B [&t|l #1 ,...l£.) B B(V 1 ,.»lft|6 J p] = d(k-k'), 

Tr B [b u \l fll ...l fj ) B B (l fl ,..A pl \bl]^f2 5 (k-Pi)S(k'-pi), (46) 

2=1 

and long time approximation, we find 

(r> + 1 )d j sin 2 tl>r!f)_t 

Ps j(t) = \n) u w <n| + ^f^|n + 1). w <n + 1| £ |/(^)| 2 j 

/m i=i (—3 — ) 

+^l»-i>. -<»-H L^^l/K-)! 22 ^^. < 47 > 

which gives the transition probability for |n) w — > |n— l) w and |n) w — >■ 
respectively as follows 

r„^„_! = Tr s [\n - l) w w (ra - l|p s /(t)] = l/M| 2 > 

m 

r_ +1 = Tr s [|n + 1> W >+ l|p sJ (t)] = (n + 1)7r ^ |/H| 2 ^^ pl -^). 



m i=i 



(45 



Specially for the choice (|22|) . we have 



r 



n— >n— 1 



m 



^'"^b^-S^-"'- (49) 

Another important case is when the reservoir has a Maxwell-Boltzman dis- 
tribution so let p/(0) = \n) u w (n\ ® p# where 
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Pb = — e v -h b , then by making use of following relations 

Tr B [b^ B bg}=Tr B [btp T B bt] = 0, 

Tnhplb^ = 5 -%^, 

eKT — 1 
5{k-k')e 



Tr B [by B b g ] = ^^>* K \ (50) 
eKT — 1 



we can obtain the density operator p s i(t) in interaction picture as 
p sI (t) : = Tr B [p/(*)] = |n) w w (n| 



|n+ l) w w<" + l| / d fc 



nw, . / ,+ °° , s , |/(u;r)| 2 e-^ sin 2 ( %~" ) t 



2m 1 '7-oo e ^_! (^)2 

which accordingly gives the following transition probabilities in verey long 
time as 



r n _> n _i = Tr s [|n - l)w u(n- l\p s l(t)] 
Tr s [|n + l) w u (n+l\p.i(t)) 



4n 2 cu 3 nt\f(uj)\ 2 e 



KT 



m eKT — 1 
An 2 u 3 (n+l)t \f(u)\ 2 



m ei<T — 1 

(52) 



substituting (j22j) in these recent relations we find 

n(5te^r 



v 

L n—>n- 



1 



m(eKT — 1 



m{eKT — \j 

So in very low temperatures the energy flows from oscillator to the reservoir 
by the rate r n ^ n _! i— > ^ and no energy flows from reservoir to oscillator. 
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4 Quantum field of reservoir 

Let us define the operators Y(x,t) and Hy(%, t) as follows 

/+oo ij^h — ' 
-=_(6^(t)e^ + 4(t)e-^), 
-oo A /2(2 7 r) 3 ^ 

/ + OO / /j_ -. 

-oo d3fc V 2(2^ ( ^ ( * )e " ^ )e ^' (54) 

then using commutation relations (j!5|) . one can show that V(x, £) and IIy(x, £), 
satisfy the equal time commutation relations 

[Y(x,t),U Y {x',t)} = i5{x-x'), (55) 

furthermore by substituting bk(t) from (EOJ) in (p2J, we obtains 

an y (f, £) 



/ + OC / qi_ 



^ w v y v ' 7-oo Y 2(2? 

BY r+°° f(u)A 

U Y (x, t) = — - 2q(t)Q(x), Q(x) = Im / d 3 k . n k) e~ ik - x , 
at J-oo ^/2(27r) 3 ^ 

(56) 

so Y(x, t) satisfies the following source included Klein-Gordon equation 

^ - V 2 Y = 2q{t)Q{x) + 2q(t)P(x), (57) 

with the corresponding Lagrangian density 

1 BY 1 _ _ BY 
£ = - -VF.VF - 2qQ(x) — + 2qP(x)Y (58) 

Therefore the reservoir is a massless Klein-Gordon field with source 2qQ(x) + 
2qP(x). The Hamiltonian density for (|57|) is as follows 

N = Y 2 ^ +-|Vlf -2qYP, (59) 

and equations (jS^jl are Heisenberg equations for Y and Ely. If we obtain 
and b? from (l5"H) in terms of 1" and Ily and substitute them in Hr defined 

k - 

in (jl4]) we obtain 

/•+oo 772 1 

Hb = J d 3 ku % blb % = ^ + -|VYf. (60) 
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5 Concluding remarks 



By generalizing Caldeira-Legget model to an environment with continuous 
degrees of freedom, for example a Klein-Gordon field, a new minimal coupling 
method introduced which can be extended and applied to a large class of dis- 
sipative systems consistently Such method applied to a quantum damped 
harmonic oscillator as a prototype of these models, with a dissipation term 
proportional to velocity. Some transition probabilities explaining the way 
energy flows between subsystems obtained. Choosing different coupling func- 
tions in (JTfij) we could investigate another classes of dissipative systems. 
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